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We report a transport study of exfoliated few monolayer crystals of topological insulator Bi2Se3 in 
an electric field effect (EFE) geometry. By doping the bulk crystals with Ca, we are able to fabricate 
devices with sufficiently low bulk carrier density to change the sign of the Hall density with the gate 
voltage V g . We find that the temperature T and magnetic field dependent transport properties in 
the vicinity of this V g can be explained by a bulk channel with activation gap of approximately 50 
meV and a relatively high mobility metallic channel that dominates at low T. The conductance 
(approximately 2 x 7e 2 //i), weak anti-localization, and metallic resistance-temperature profile of the 
latter lead us to identify it with the protected surface state. The relative smallness of the observed 
gap implies limitations for EFE topological insulator devices at room temperature. 

PACS numbers: 73.20.-r, 73.20.Fz, 73.23.-b, 73.63.-b 



Topological insulators (TIs) are a new phase of mat- 
ter which is electrically insulating in the bulk but has 
unusual conducting surface states (SS) [iHEj]- The SS 
are spin polarized, protected from scattering by non- 
magnetic impurities, and have an approximately linear 
energy-momentum dispersion at low energy. These prop- 
erties make them relevant for applications ranging from 
improved spintronic devices to producing analogs of ex- 
otic high energy particles potentially useful for quantum 
computing [6]. Particular interest has been focused on 
3D TIs based on Bi, where the topologically non-trivial 
surface 2D electron gas has been positively identified by 
angle-resolved photoemission spectroscopy (ARPES) in 
Bi-Sb 0, Bi 2 Se 3 [8], and Bi 2 Te 3 @. Proposals for real- 
izing devices that harness the aforementioned properties 
of the SS most often involve electrical transport [Iol4l4). 
Thus far, characterizing or even identifying the SS in such 
experiments has proven difficult. Because transport is a 
bulk sensitive measurement, even a small conductivity 
from imperfections in the bulk overwhelms the surface 
contribution because of the geometric advantage of the 
former. The materials obstactles to placing the chemical 
potential \i in the bulk band gap of the TI are at present 
the most serious obstacle to realizing TI devices. 

Recent measurements in high magnetic fields have 
identified the surface state as a parallel conductance 
channel with the bulk [l5l-[l7j . Several experiments have 
been recently undertaken to circumvent the bulk contri- 
bution to transport in TIs by preparing Bi2Se3 in thin 
film form [18l-[21j. However, in both bulk crystals and 
nanocrystals of Bi2Se3, \i is usually pinned to the con- 
duction band (CB) due to the presence of Se vacancies 
pH |22| . This large remnant electron density is too large 
to remove by the electrostatic gating techniques available 
for thin films. Thus, while thin film experiments have 
been able to show signs of surface phenomena, namely 
remnants of weak anti- localization [18|, [2l| , basic trans- 
port properties of TIs have remained elusive as the bulk 
carriers cannot be effectively removed. 
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FIG. 1: (a) R vs. V g at selected T between 10 and 120 K for 
sample SI. The inset is a photo of device SI with scale-bar (4 
H- (b) Cuts of R vs. T with V g fixed at -100, -50 and +25 
V. Gating results for samples S3, S4, and S5 are shown in 
panels (c), (d), and (e). An optical image of a cleaved crystal 
with various step heights with thickness labels measured by 
atomic force microscopy (in nm) is shown in the inset of (e). 



Here we report experiments on Ca passivated crystals 
of Bi2Se3 in a field effect transistor device that allow us to 
successfully suppress the bulk states and investigate the 
transport properties of the SS. We mechanically exfoli- 
ated crystals of Ca doped Bi2Se3 (up to 0.5% substituted 
for Bi) onto doped silicon wafers coated with 300 nm 
SiC>2. The Ca doping results in hole doped bulk crystals 
[23l . l24j , which after fabrication of contacts via e-beam 
lithography (Cr/Au) in the Hall geometry have \i rela- 
tively close to the bulk band gap on the CB side. In sam- 
ples with uniform thickness a negative V g larger than 
-150 V can be applied without triggering breakdown. For 
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d < 10 nm, from electrostatic arguments this suffices to 
move \i deep into the gap for bulk carrier densities 10 18 
- 10 19 cm" 3 . 

The use of 300 nm Si02 allows optical identification 
of thin crystals. It was found that for crystals thinner 
than ~ 30 nm the observed color of the cleaved crystals 
darkens toward the color of the substrate. This is sim- 
ilar to what has been found in graphene and other 2D 
crystals [25j. Cross correlating the optical images with 
atomic force microscope (AFM) images, we were able to 
generate a mapping between color and thickness useful 
for identifying thin crystals optically within an approx- 
imate thickness ± 2 nm. An optical image of a crystal 
with varying thickness measured by both AFM and opti- 
cal microscopy is shown in the inset of FigQJ. We study 
crystals down to 5 nm; ARPES experiments on thin films 
of Bi2Se3 have shown the SS persist to this thickness (26^ . 

Table U shows parameters for five different devices. De- 
vices of thickness t < 20 nm exhibit similar characteris- 
tics. For these devices an applied gate voltage can change 
the Hall response from electron-like to hole-like at the 
neutral voltage Vn- Each exhibits a maximum at low T 
in the resistance per square R mdiX of 2 — 3 kfi, as shown 
for samples SI and S3 in Fig. QJi and QJ, respectively. 
For thicker crystals such as devices S4 and S5, we reach 
the threshold for breakdown of the Si02 gate dielectric 
(typically -180 V to -200 V) without observing such a 
peak (see Fig. [T]i and The progressively thicker 

crystals show systematically less response to the applied 
gate voltage. It is found that the act of cleaving and de- 
vice fabrication reduces the carrier mobility and increases 
the electron density. Device S5 (see Fig. le and Table 
I) has a 0.3 K mobility approximately 1090 cm 2 V _1 s _1 
with 6 x 10 18 e~ cm -3 was cleaved from a bulk crystal 
measured to have a mobility 7000 cm 2 V _1 s _1 with 2 
x 10 18 e~ cm -3 . It is unclear if the cleaving process or 
lithographic process is responsible for this degradation. 
From here we focus on the results of samples SI and S2. 
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10 ± 2 
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30 ± 5 
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TABLE I: Sample parameters, t is the sample thickness esti- 
mated from an optical image of the cleaved crystal (see inset 
of Fig [T^) . Vn is the gate voltage required to reach zero Hall 
voltage at 5 K. fib is the Hall mobility in the conduction band, 
calculated at zero gate voltage. i? ma x is the maximum resis- 
tance recorded over the measured gating range. 



Figure [TJi displays traces of the zero-field resistance 
R = R xx vs. V g at selected T from 10 to 120 K (sam- 
ple SI). For V g > -55 V, R is nearly T independent, 
characteristic of a bad metal in which disorder scattering 
is dominant (Fig. [TJd). Surprisingly, if fi is tuned lower 
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FIG. 2: (a) Hall resistance R yx vs. 
K (bold curve) to 130 K for sample SI. The sketch (inset) 
shows the strong band-bending induced by V g <-80 V. The 
thermal excitation of carriers is indicated by red arrows (f(E) 
is the Fermi-Dirac distribution), (b) Fit of the conductance G 
at V g =-im V (solid circles) to Eq. [IJ with G s (0) = 13.9 e 2 /h. 
Values of n&(T) (open circles) inferred from R yx and AV g are 
plotted as open circles, along with the fit to the activated 
form n b = n bo e" A/T (A = 640 K). 



(V g < -55 V), R displays a power-law variation R ~ T a 
(a=2.8) that extends to 110 K, where R displays a broad 
maximum. We show below (Eq. [T]) that both the gate- 
induced transformation (from a T-independent R to one 
with power-law variation) and the existence of the broad 
peak in R are accounted for by a 2-band model com- 
prised of a bulk band in parallel with a surface channel 
with higher mobility at this large negative gate voltage. 
When V g is near -55 V, fi appears to be at the CB band 
edge. 

The Hall resistance R yx confirms that the carriers are 
n-type for V g > -55 V (Fig. [2^). R yx is T independent, 
consistent with a bulk metal. Additionally, R yx (V g ) is 
linear as we would expect from a simple electron band be- 
ing emptied by the gate (induced charge rii n d = CV g /e, 
where C the capacitance per unit area of the gate dielec- 



tric). As V g is lowered below -55 V, the curves of Ryx vs. 
V g display a strong dependence on T. Below 30 K, R yx 
changes sign near V g = -80 V. We can understand this 
in terms of fi moving below the SS Dirac point or near 
the VB edge. It is important to note that the gating does 
not rigidly shift \i in the device. For samples with d <C D 
(the depletion length), the bands are strongly bent up- 
wards when V g is large and negative (inset in Fig. [2^). 
As ji is grounded at the drain, it falls within the gap in- 
side the crystal. While the specific shape of the bending 
is difficult to calculate, we can adopt the most simple 
limit by treating the system as a metal - (doped) semi- 
conductor interface in the full depletion approximation. 
The Poisson equation for the potential </>, charge density 
n (estimated as 10 19 e~ cm -3 from R yx ), and dielectric 
constant e as a function of the distance from the interface 
z is then </>(z) = (en/2e)(D 2 — (D — z) 2 )). The curva- 
ture of the band is that depicted in Fig. [2^ for a 10 nm 
crystal that should be preserved in more realistic calcu- 
lations, which will be important for analyzing the role of 
thermally activated carriers. 
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We assume that, when fi is inside the gap, the observed 
conductance G = 1/R is the parallel combination 

G(T) = G s (T) + G b (T), (1) 

where G S (T) = n s eji s is the surface conductance, with 
n s the two-dimensional (2D) density and ji s — T _Q! 
the mobility (e is the electron charge and a = 2.8). 
To account for disorder scattering, we write G S (T) = 
G s (0)/[1 + AT a ], and treat the zero-T surface conduc- 
tance G s (0) and A as adjustable parameters. The sec- 
ond term G b (T) = ntefib is a bulk conductance result- 
ing from thermal excitation of carriers into the CB de- 
scribed by the density ra&(T) = nboe~ A / kBT (arrows in 
inset, Fig. [2Jl). The mobility /i& at the bottom of the 
CB is taken to be a constant (disorder-scattering domi- 
nant). We note that, in addition to \i originally being in 
the CB, the curvature of the band bending in the case 
of a remnant electron pocket favors excitation in to the 
CB for devices of finite d. Integrating the excited car- 
riers n exc ~ J d e(e- A ^/^ T - e~ A ^/ k - T )dz (A CB ,vb 
is the energy difference between the VB and CB band 
edges and n) , we can expect the electron-like carriers to 
dominate due to the direction of the band bending. 

Equation [T] gives a good fit to G vs. T (Fig. [2]o) . From 
the fit, we obtain G s (0) = 13.9 e 2 /h (A = 2.5 x KT 6 ). 
This implies a conductance of approximately 7e 2 /h per 
surface. As T increases from 10 K, inelastic scattering 
strongly reduces G S (T). Meanwhile, excitation across 
the gap A = 530±30 K leads to a rapid growth in the 
bulk G b . This reverses the decreasing trend of G(T) at 
110 K, as shown in Fig. At 130 K, the 2 terms, G b (T) 
= 3.65 e 2 /h and G S (T) = 4.5 e 2 /h, become roughly com- 
parable. Thus, we show that despite the large intrinsic 
gap in Bi2Se3 (300 meV), we must go below 100 K for 
the device to be dominated by surface conduction. One 
likely origin of this reduction in energy gap is the band 
bending depicted in Fig. [2^, which efflctively narrows 
the gap upon increased gating. This implies limitations 
for accessing the SS at room temperature. 

To measure A more accurately, we turn to the Hall 
results in Fig. [2^. We may obtain n^T) by studying 
how R yx depends on both V g and T when \i is in the 
gap. At the lowest T (bold curve at 10 K), R yx is at 
the zero-crossing when V g = -80 V. As T is raised with 
V g fixed, the balance condition is unaffected until T >80 
K, whereupon R yx swings to negative values, with the 
thermal activation of carriers into the CB. However, the 
balance point can be restored if V g is made more negative 
to draw more holes into the crystal (e.g., V g —> -98 V 
at 100 K). From the gate-shift AV g needed at each T, 
we calculate the activated density as n^(T) = AV g C/e. 
A fit to n b (T) yields a gap of 640 K (Fig. Eb). This 
excitation supports the band bending proposed in Fig. 
2a. Combining the value — 6 x 10 12 cm -2 at 130 K 
with the value of G b (T) inferred from Eq. [TJ we obtain 
the mobility /i& — 147 cm 2 /Vs (at the bottom of the CB). 

Hence, the gate-driven transformation to a conductor 
with strong T dependence in both R and R yx is explained 
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FIG. 3: (a) Conductance G(T, 0) in zero H (upper trace) 
with G(T,H) at 14 T (lower) at T = 5 K for sample SI. 
In the H = trace (superposition of 3 traces), retraceable 
fluctuations are resolved, (b) Amplitude of the conductance 
fluctuation 5G — G — (G) (with (G) a smooth background) 
(c) The T dependence of the rms value 5G r ms (solid circles 
in inset) fits well to T -0-5 (dashed curve). 



quantitatively as the opening of a small gap A —50 meV 
between \i and the CB. Depletion of the CB carriers ex- 
poses a high- mobility channel G s . The possibility that 
this channel is a bulk impurity band can be excluded 
because electrons therein generally have even lower mo- 
bilities (1-10 cm 2 /Vs) than obtained for 

The enhanced mobility for G s is confirmed by the 
transverse magnetoresistance (MR). Figj3k compares the 
zero-field conductance G(T, 0) at 5 K with the conduc- 
tance G(T, H) measured in with H||c at 14 T as a func- 
tion of V g (we suppress T in G(T, H) hereon). The semi- 
classical MR expression G(H) = G(0)/[1 + (/w#) 2 ] 
may be used to determine the average mobility \i ave . At 
V g -20 V (/i in the CB), the reduction in G by 25% 
gives fi ave —380 cm 2 /Vs, consistent with the Hall analy- 
sis (see Table I). The larger reduction (60%) in the gap 
region (V g —-80 V) gives \i ave = 850 cm 2 /Vs, which we 
identify with \i s . From G s , we can estimate n s — 4, 10 x 
10 12 cm" 2 at V g = -80,0 V. 

At low T, fluctuations in G are observed if H or V g 
is swept. Fig. [3)3 displays the trace of fluctuating com- 
ponent SG = G — (G) vs. V g , where (G) is a smooth 
background. The amplitude of SG is nominally uniform 
over the entire range —100 < V g < 20 V. The ampli- 
tude (Fig. [3b) follows a T -1 / 2 scaling consistent with a 
2D electron system with a dephasing length smaller than 
the system size [27|. Recent observation of reproducible 
fine structure driven by H in both transport and STM 
measurements have been reported 0, [28| . Here, for Vg 
= (Fig. 2^), we observe in Sample S2 (R(V g ) shown in 
Fig. [4jd) a sharp, narrow peak at H = flanked by mesa- 
like structures in addition to the fluctuations. However, 
the large amplitude fluctuations (rms amplitude up to 
6e 2 /h) that we found in bulk samples of the same ma- 
terial previously [24] do not survive in these exfoliated 
devices. The reason for this and the origin of the de- 
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FIG. 4: Panels (a) and (b) show conductance fluctuations for 
sample S2 in G vs. H at 0.3 K with V g fixed at V and at -90 
V. (c) Low- if magnetoconductance at V g = V and V g =-100 
V. The solid lines are fits to Eq. [2] (d) Fit parameter A shows 
a sharp maximum near the charge- neutral point. Dashed line 
shows the value A = 1/tt predicted for dominant spin-orbit 
coupling (see text). R(V g ) at 0.3 K is also shown. 



tailed structures observed (including the suppression by 
H shown in Fig. [3h) deserve further study. 

The sharp anomaly in G at H = in Figs. 2^ and0)3 
is consistent with the quantum correction arising from 
anti-localization given by [29] 



AG XX (H) = A- 



i #° / f 1 j. H » 



(2) 



where Hq is the dephasing magnetic field, i/j is the 
digamma function, and the parameter A is positive for 
antilocalization. If spin-orbit coupling is strong, theory 



predicts A = l/2n [29j, or 2 A = 1/tt if we have 2 sur- 
faces, as here. Confining the fit to fields below 0.4 T, 
we show 2 traces in Fig. [4fc at V g = and -100 V. Fits 
across the full range of V g are shown in Figure [4jT. At 
Vat, A peaks at 0.38, which is within 20% of the expected 
value 1/tt. Reports of similar experiments on epitaxially 
grown films report approximately half of this value, con- 
sistent with one surface state [18|, [21( . We hypothesize 
the disparity may arise from the different nature of the 
strained interface of the epitaxial films and the exfoliated 
crystal here. We also find a peak in the dephasing field 
Ho close to Vn indicating a reduction in the dephasing 
length from 220 nm at zero V g to 67 nm at -90 V. 
The value at zero gate is consistent with those reported 
in epitaxial films of similar thickness [2l|, while the re- 
duction in has been suggested as evidence for reduced 
screening and increased electron-electron interaction ef- 
fects in the low density regime [l8j. The observation of 
A ~ 1 /tt along with the sign change in R yx , supports our 
main conclusion that, at Vat, [i is well inside the bulk 
band gap. Notably, even with \i in the CB, a signature 
of weak anti-localization persists. 

Electrostatic control of \i provides a powerful means for 
comparing the bulk and surface conductances in chem- 
ically doped TIs. We anticipate that replacing the Au 
electrodes with magnetic or superconducting materials 
will enable us to probe exotic states that have been pro- 
posed |lCH14j . The current devices allow access to surface 
dominated transport for these experiments below 100 K. 
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